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ABSTRACT 

We have developed a Monte Carlo technique to test models for the true power spectra 
of intermittently sampled lightcurves against the noisy, observed power spectra, and 
produce a reliable estimate of the goodness of fit of the given model. We apply this 
technique to constrain the broadband power spectra of a sample of four Seyfert galaxies 
monitored by the Rossi X-ray Timing Explorer (RXTE) over three years. We show that 
the power spectra of three of the AGN in our sample (MCG-6-30-15, NGC 5506 and 
NGC 3516) flatten significantly towards low frequencies, while the power spectrum of 
NGC 5548 shows no evidence of flattening. We fit two models for the flattening, a 'knee' 
model, analogous to the low-frequency break seen in the power spectra of BHXRBs in 
the low state (where the power-spectral slope flattens to a = 0) and a 'high-frequency 
break' model (where the power-spectral slope flattens to a = 1), analogous to the high- 
frequency break seen in the high and low-state power spectra of the classic BHXRB 
Cyg X-1. Both models provide good fits to the power spectra of all four AGN. For 
both models, the characteristic frequency for flattening is significantly higher in MCG- 
6-30-15 than in NGC 3516 (by factor ~ 10) although both sources have similar X-ray 
luminosities, suggesting that MCG-6-30-15 has a lower black hole mass and is accreting 
at a higher rate than NGC 3516. Assuming linear scaling of characteristic frequencies 
with black hole mass, the high accretion rate implied for MCG-6-30-15 favours the 
high-frequency break model for this source and further suggests that MCG-6-30-15 
and possibly NGC 5506, may be analogues of Cyg X-1 in the high state. Comparison 
of our model fits with naive fits, where the model is fitted directly to the observed 
power spectra (with errors estimated from the data), shows that Monte Carlo fitting 
is essential for reliably constraining the broadband power spectra of AGN lightcurves 
obtained to date. 
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ical 



1 INTRODUCTION 

The strong and rapid X-ray variability observed in many 
Seyfert galaxies on time-scales of a day or less provides 
strong evidence that the X-rays are emitted close to the 
central black hole. Early efforts to characterise the X-ray 
variability, using data from EXOSAT, showed that it has 
a scale-invariant, red- noise form on time-scales from a few 
hundred seconds up to the few-days duration of the observa- 
tions (M'^Hardy & Czerny 1987, Lawrence et al. 1987) . Later 
studies of the X-ray variability properties of large samples 
of radio-quiet AGN showed that the variability amplitude 
scales inversely with luminosity (Green, M'^ Hardy & Lehto 
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1993; Lawrence & Papadakis 1993; Nandra et al. 1997). One 
possible explanation of this result is that the higher luminos- 
ity AGN contain more massive black holes and the variabil- 
ity time-scales in AGN scale with black hole mass. Intrigu- 
ingly, black hole X-ray binary systems (BHXRBs) also show 
red-noise type variability of a similar amplitude to AGN, 
on time-scales less than seconds. The similarity in X-ray 
variability properties of AGN and BHXRBs raises the pos- 
sibility that the processes causing variability in AGN and 
BHXRBs are the same and that any characteristic variabil- 
ity time-scales scale with the central black hole mass. This 
possibility can be tested by comparing the detailed X-ray 
timing properties of BHXRBs and AGN. 
Timing studies of BHXRBs are usually carried out in the 
frequency domain using the power spectrum, which shows 
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the contribution of variations on different time-scales (cor- 
responding to power-spectral frequencies) to the total vari- 
ability of the lightcurve. The power spectra of BHXRBs 
are dominated by a broadband noise component (van der 
Klis 1995). On short time-scales, the variability is charac- 
terised as scale-invariant 'red noise', producing a power-law 
power spectrum (power P(i') at frequency v is given by 
P(z/) oc where a is the power-spectral slope) of slope 
a ~ 1-2.5 (van der Klis 1995). In the 'low' state, charac- 
terised by a relatively hard X-ray spectrum (similar to that 
of AGN), the power spectrum flattens towards lower fre- 
quencies so that on long time-scales the X-ray lightcurve 
becomes 'white noise', with corresponding slope a ~ 0. For 
example, in the classic BliXRB system Cyg X-1, the power- 
spectral flattening is well described by a power-law with two 
breaJis, a high-frequency break which varies between 1 and 
6 Hz, above which the power-spectral slope varies between 
Q ~ 1.5-2.4 and a low- frequency break which varies between 
0.04 and 0.4 Hz, above which the power-spectral slope a ~ 1 
and below which the slope a ~ (Belloni & Hasinger 1990). 
In contrast, the power spectrum of the 'high' (soft energy 
spectrum) state seen in some BHXRBs (inluding Cyg X-1) 
does not flatten to zero slope; instead, the slope q ~ 1 below 
the high-frequency break extends to < lO^'^ Hz (e.g. Cui et 
al. 1997). 

In order to test the hypothesis that the X-ray variability 
of AGN is similar to that of BHXRBs over a broad range 

of time-scales, we must search for low-frequency flattening 
in the broadband power spectra of AGN. By fitting models 
with power-spectral breaJis to the AGN power spectra and 
comparing the estimated break frequencies with what we ex- 
pect if they correspond to similar breaks in BHXRB power 
spectra, we can test the possibility that the power-spectral 
shape is really the same and scales simply with black hole 
mass. If so, we expect break frequencies in AGN to be found 
at frequencies of ~ 10~^ Hz or lower, so that monitoring ob- 
servations on time-scales of weeks or longer are necessary to 
detect any flattening in the power spectrum. 
Early attempts to measure broadband power spectra of 
AGN were hampered by the sparseness of long-term archival 
lightcurves, which had to be constructed from data obtained 
by several missions (IVf^Hardy 1988). Nonetheless, some ev- 
idence for power-spectral flattening was found, but models 
for the form of the flattening could not be constrained (Pa- 
padakis & M-'Hardy 1995). 

Ideally, broadband power spectra should be measured from 
lightcurves obtained with frequent and regular sampling over 
a long duration, which previous missions were not optimised 
to do. The Rossi X-ray Timing Explorer (RXTE), which 
has just such a capability, was launched in December 1995. 
RXTE carries a large-area proportional counter array (the 
PCA) which can detect many AGN with good signal-to- 
noise in less than 1000 s, but most importantly, RXTE can 
slew rapidly so that it may monitor many targets with fre- 
quent 1 ks snapshots. 

We have monitored a sample of 4 Seyfert galaxies (MCG-6- 
30-15, NGC 4051, 5506 and 5548) with RXTE since 1996, 
in order to measure their broadband power spectra. These 

objects are known to be significantly X-ray variable and 
cover a broad range of X-ray luminosity (NGC 4051~ 
5 X 10"^ erg s"\ MCG-6-30-15 and NGC 5506~ 1.5 x 
10*^ erg s-\ NGC 5548- 5 x 10"*^ erg s"^) and presum- 



ably, a broad range of black hole masses. We describe the 
power spectrum of NGC 4051, which shows urmsual non- 
stationarity in its lightcurve (Uttley et al. 1999) in a sep- 
arate paper (Papadakis, M'^Hardy &: Uttley, in prep.). In 
this paper, we present a power-spectral study of the re- 
maining three objects in our sample using data from RXTE 
cycles 1, 2 and 3, also including the excellent lightcurves 
obtained as part of a separate power-spectral study of the 
Seyfert 1 galaxy NGC 3516 (luminosity- 1.5 x 10"*^ erg s"^), 
by Edelson & Nandra (1999). We describe data reduction 
and present the lightcurves in Section 2. 
The estimation of the underlying power-spectral shape from 
lightcurves which are discretely (and possibly unevenly) 
sampled is hampered by the distorting effects of aliasing 
and red-noise leak;. A further serious problem is that the 
measured power spectra are intrinsically noisy, and reliable 
errors on the power in each frequency bin cannot be es- 
timated from the data (especially at low frequencies), due 
to the small number (< 20) of power-spectral measurements 
made in each frequency bin. In Section 3, after presenting the 
observed power spectra, we describe these problems, which 
previous efforts to constrain the shape of the broadband 
power spectrum of AGN using RXTE data (e.g. Edelson & 
Nandra 1999, Nowak & Chiang 2000) have not accounted 
for. 

To overcome the difficulties in estimating the true power- 
spectral shape, we have developed a method which we call 
PSRESP, based on the response method (Done et al. 1992) 
which uses Monte Carlo simulations of lightcurves to take 
account of the distorting effects of sampling and to estimate 
uncertainties, allowing us to test various power-spectral 
models against the data. We describe psresp in Section 4, 
and apply it to the lightcurves of our sample of Seyfert 
gala^xies in Section 5, in order to test for flattening in their 
broadband power spectra and constrain simple models for 
describing any flattening we see. In Section 6, we compare 
our results with those obtained by naively fitting the ob- 
served power spectrum (without taking proper account of 
errors and the distortion due to sampling), use our power 
spectral measurements to estimate the black hole masses of 
the AGN in our sample and discuss some of the implica- 
tions of our results, before making concluding remarks in 
Section 7. 

2 THE LIGHTCURVES 

2.1 Observations and data reduction 

We use data obtained with the PCA on board RXTE cov- 
ering a three year period during observing cycles 1-3 when 
the PCA gain setting was constant, so that count rate m&ar 
surements provide a simple measure of observed flux. Of the 
three instruments on board RXTE, only the PCA is sensitive 
enough to allow us to make an accurate flux measurement 
for our targets with the 1 ks snapshots which our monitoring 
consists of. The PCA consists of 5 Xenon-filled Proportional 
Counter Units (PCUs) , numbered to 4 which are sensitive 
in the 2 60 keV energy range and contribute to a total effec- 
tive area of 6500 cm^. Since launch, discharge problems have 
meant that one or both of PCUs 3 and 4 are often switched 
off (this problem extended to PCU 1 in March 1999 but we 
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Table 1. Lightcurve details 



MCG-6-30-15 
NGC 5506 
NGC 5548 
NGC 3516 



Long-term 
Start Stop 
8 May 1996 2 Feb 1999 
23 Apr 1996 2 Feb 1999 



Intensive 
Start Stop 
23 Aug 1996 29 Sep 1996 
8 Aug 1996 19 Sep 1996 



Long-look 

Start Stop Exp 

03:31 UT 4 Aug 1997 12:34 UT 12 Aug 1997 332 

04:45 UT 20 Jun 1997 12:33 UT 9 Jul 1997 93 



23 Apr 1996 22 Dec 1998 26 Jun 1996 8 Aug 1996 12:41 UT 19 Jun 1998 07:01 UT 24 Jun 1998 99 
16 Mar 1997 28 Dec 1998 16 Mar 1997 30 Jul 1997 00:14 UT 22 May 1997 05:37 UT 26 May 1997 249 



The table shows the start and stop times of the lighteurves used in this work (except 2nd NGC 3516 long-look - 
Also given is the useful exposure time in ks (Exp) for each long-look observation. 



see text for details. 
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Figure 1. 2-10 keV monitoring lighteurves of MCG-6-30-15, NGC 5506, NGC 5548 and NGC 3516. CR units are count s'^ . 



do not include this later data here). Despite the loss of up 
to two PCUs during our observations, we are easily able to 
obtain sufficient signal-to-noise in a single snapshot for our 
purposes {S/N > 20a). 

In order to efficiently measure a power spectrum over the 
broadest range of time-scales while minimising the neces- 
sary observing time, we monitored our targets using several 
different schemes, each designed to measure the power spec- 
trum over a different frequency range. In 1996, we observed 
MCG-6-30-15, NGC 5506 and NGC 5548 twice daily for 
~ 2 weeks followed by daily observations for ~ 4 weeks and 
then weekly for the remainder of the year. During the fol- 
lowing two years, we observed our targets every two weeks. 
NGC 3516 was monitored as part of a separate study with 
broadly similar goals to our own (Edelson & Nandra 1999) . 
Here we use public archival data from this campaign, in- 
cluding an intensive period of monitoring every 12.8 h for 
^4 months duration, and long-term monitoring at 4.3 d in- 
tervals from March 1997 until the end of 1998. The start 
and end dates of all the lighteurves are shown in Table ^ 
We measure variability on short time-scales using 'long- 
look' observations, quasi-continuous observations of dura- 



tion ~days, which we obtained ourselves (NGC 5506) or 
from the RXTE public archive (MCG-6-30-15, NGC 5548 
and NGC 3516). The details of these observations are also 
summarised in Table |l]. Not shown in Table ^ are details of 
a second long-look observation of NGC 3516, obtained from 
08:00 UT 13 April to 16:13 UT 16 April 1998 (148 ks useful 
exposure), which we also include to maximise the definition 
of the power spectrum of NGC 3516 at high frequencies. 
Unfortunately, the NGC 5506 long-look is too sparsely sam- 
pled (spread over a 20 day period) to be useful for mea- 
suring the power spectrum except at the highest frequencies 
(~ 10~^ Hz). Therefore, in order to measure the power spec- 
trum of NGC 5506 in the lO'^'-lO"^ Hz range, we use an 
archival EXOS AT ME lightcurve, of ~ 230 ks continuous du- 
ration, obtained during 24-27 January 1986 and originally 
described by M'^ Hardy & Czerny (1987). The energy range 
sampled by EXOS AT (1-9 keV) is comparable to the 2- 
10 keV range which we will measure with RXTE, so that 
the normalised power spectrum should have a similar shape 
and amplitude to that measured by RXTE, if the power- 
spectral shape is stationary on time-scales of a decade (see 
Section O). 



© 0000 RAS, MNRAS 000, 000-000 



4 P. Uttley, I. M. Hardy and I. E. Papadakis 



30 - MCG-6-30-15 

K 20 - . ■ . • • . . • . . - 

" :. . .' .■*.-••.•• • • . . 

10 ■ ■ ■ ■ . • , , • • : 

' I I I I I I I I I I I I I I I I I I I I I I I I I i I I I I i I I I I I I I I I I ' 
NGC 5506 

40 - 

K . . ■ . . . 

o . ' ■ ' . . ■ ■ . ■ . 

20 - ■ . ■ ' ■ ■ ■ • ' 

- I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I - 
30 r NGC 5548 t 

K 20 . ....... ■ ■ - 

u ■ ■ • • . ^ _ ^ . . 

10 - 

: I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I : 

30 : NGC 3516 : 

K 20 : . ■ • I ■ ' I ' . . : 

10 : ' ■ .J 

L- 1 I I I I I I I I \ I I I I \ I I I I I I I I I I I I I I \ I I I I \ I I I I I 1— 

5 10 15 20 25 30 35 40 

Time (days) 

Figure 3. 2-10 keV intensive monitoring lightcurves of MCG-6-30-15, NGC 5506, NGC 5548 and NGC 3516. 
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Figure 4. 2-10 lieV long-looli lightcurves of MCG-6-30-15, NGC 5506 {EXOSAT 1-9 IceV), NGC 5548 and NGC 3516. 



We reduce all RXTE data using ftools v4.2. Because 
PCUs 3 and 4 are often switched off, we only use data 
from PCUs 0, 1 and 2. We extract data from the top layer 
only (to minimise background relative to source counts) and 
make lightcurves in the 2-10 keV channel range correspond- 



ing to absolute channels 7-28. Wo exclude data obtained 
within and up to 20 minutes after SAA maximum and data 
obtained with earth elevation < 10°, target offset > 0.02° 
and electron contamination > 10%. We estimate background 
lightcurves using the L7 background model. 
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Figure 2. Comparison of ASM and PCA lightcurves of 
NGC 5548. The lightcurves are averaged in 28 d bins. ASM er- 
ror bars represent standard errors. No errors are determined for 
the PCA Hghtcurve, although the small spread of points on short 
time-scales in the original lightcurve (Fig. |^) suggests that the 
uncertainties are small. 



We show the long-term monitoring lightcurves in Figure |l| 
The annual gaps lasting ^ 6-8 weeks in the MCG-6- 
30-15 and NGC 5506 lightcurves correspond to periods 
when sun-angle constraints prevent RXTE from pointing 
at these objects. Strong variability can be seen in all four 
Hghtcurves, and long-term trends are particularly apparent 
in the lightcurves of NGC 5548 and NGC 3516. 
It is important to note that the quality of these long-term 
monitoring lightcurves is far superior to that obtainable with 
the All-sky monitor (ASM) on board RXTE which, although 
excellent for monitoring bright sources, is subject to large 
systematic errors when observing faint sources like the AGN 
we study here. This is apparent if we compare the 28-day 
averaged ASM and PCA lightcurves of NGC 5548 obtained 
over the same period (see Fig. ^). The ASM lightcurve looks 
very different to the PCA lightcurve, therefore ASM data 
should not be used to measure the low-frequency power spec- 
tra of faint sources (< 0.5 ASM count s~^). 

A close-up look at the period of intensive (twice-daily 
and daily) PCA monitoring can be seen in Figure |^, which 
is plotted in terms of days since the start of each inten- 
sive monitoring period. The NGC 3516 intensive monitoring 
lightcurve is cut short so that the lightcurves are of simi- 
lar length for comparison purposes (see Edelson & Nandra 
1999, for the full lightcurve). Occasional short gaps in the 
lightcurves are due to observations excluded due to our data 
extraction criteria. Significant variability on time-scales of 
days can be seen in all four lightcurves, but MCG-6-30-15 
shows the strongest variations on the shortest time-scales. 
We show the long-look lightcurves in Figure |^, binned to 
512 s resolution. For comparison purposes, we plot simi- 
lar lengths of lightcurves and cut off more than half of the 
MCG-6-30-15 lightcurve (which can be seen in full in Lee 
et al. 1999) . We plot only the April 1998 long- look observa- 
tion of NGC 3516 (see Edelson & Nandra 1999 for the ear- 
lier long-look lightcurve). We show here the continuous iJJ^- 
OSATlightcurve of NGC 5506 for comparison purposes (see 
Lamer, Uttley & M'^Hardy 2000 for the RXTE lightcurve). 



On short time-scales, it can be seen that the MCG-6-30- 
15 and NGC 5506 lightcurves look similar and show quite 
strong, rapid variability. On the other hand, NGC 5548 and 
NGC 3516 show slower, more gradual trends. 



2.2 Background subtraction and source 
contamination 

Because the PCA is not an imaging instrument, the con- 
tribution of background to the lightcurves must be mod- 
elled. Discrepancies between the model background and the 
real background might then contaminate the background- 
subtracted lightcurves. A further source of contamination 
may be due to other, reasonably bright sources in the field 
of view. We now briefly consider the possible contribution 
of this contamination to our lightcurves. 
Edelson & Nandra (1999) use offset pointings to show that 
the average discrepancy between the L7 background model 
and the measured background in the 2-10 keV band is sig- 
nificant (0.87 count s~^) but varies little (noise subtracted 
RMS 0.39 count s~^) and so introduces little power into 
the measured power spectrum. Moreover, the variations in 
this background error occur only on long time-scales (weeks) 
where the source variability is stronger, so we do not ex- 
pect any spurious power introduced by these variations 
to be significant compared to the power intrinsic to the 
source. However, Uttley et al. (1999) show that spectra of 
NGC 4051, obtained simultaneously while it was very faint 
(1.3x 10"^^ erg cm"^ s"^) by RXTE and the imaging MEGS 
intruments on board BeppoSAX are in good agreement with 
one another, implying little background offset (since the to- 
tal 2-10 keV source count rate observed by RXTE in 3 PCUs 
was only 0.4 count s~^). This discrepancy with the signifi- 
cant background offset observed in NGC 3516 implies that 
the offset is dependent on the source being observed, and 
hence may be associated with spatial fluctuations in the cos- 
mic X-ray background or other faint sources in the RXTE 
fleld of view. Any small constant offset due to inaccurate 
background modelling will only affect the normalisation of 
the power spectrum by a relatively small amount (once it 
has been normalised by squared mean flux, see Section 3.1) 
and will not affect the shape of the power spectrum at all, 
so we do not consider it further in the cases of MCG-6-30- 
15, NGC 5506 and NGC 3516. We note however that the 
observations of NGC 5548 suffer a minor complication, in 
that the field of view also contains the bright BL Lac ob- 
ject IE 1415. 6-f 2557, offset 0.5° from NGC 5548. Chiang 
et al. (2000) conducted separate pointings at this source 
and found that its contaminating contribution to the mea- 
sured 2-10 keV PCA count rate (for 3 PCUs) of NGC 5548, 
after allowing for the effects of the PCA collimators, was 



only ~ 2 count s~ (about 10% of the total measured 
count rate). The contaminating flux was estimated to vary 
by < 0.8 count s~^ in two months so that, assuming that 
there is not much stronger variability on longer time-scales, 
IE 1415.6-1-2557 should not contribute significantly to the 
low-frequency power measured from the RXTE lightcurve. 
However, in order to take account of the contaminating con- 
tribution to the mean flux level of the NGC 5548 lightcurve, 
we shall subtract 2 count s~^ from the measured 2-10 keV 
mean flux level of NGC 5548 for the purposes of power- 
spectral normalisation. 
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3 THE POWER SPECTRA 

3.1 Measuring the raw power spectra 

To obtain the power spectrum of a discretely and possibly 
unevenly sampled light curve f{ti), of length A'' data points, 
we first subtract the mean flux fi from the lightcurve (to 
remove zero-frequency power) and then calculate the modu- 
lus squared of its discrete Fourier transform at each sampled 
frequency v (e.g. Deeming 1975); 



\Fn{i^) 



^f{ti) cos(27ri/ti) 1 + i'^fiU) sm{2Tviyti) 



Note that the frequencies sampled by the discrete 
Fourier transform occur at evenly spaced intervals, 
i/niin, 2fmin, Sumin, ....i^Nyq, where Umin is equal to (where 
T is the total duration of the lightcurve, i.e. T = In — ti) 
and the Nyquist frequency !^Nyq = {2T/N)~^ . We obtain 
the power P{i') by applying a suitable normalisation to 
\Fn{i')\'^- Throughout this work we apply the fractional 
RMS squared normalisation, 

which is commonly used in measuring XRB power spectra 
and has the desirable property that integrating the power 
spectrum over a given frequency range, Vi to 1^2 yields the 
contribution to the fractional RMS squared variability (i.e. 
/fj.^) of the lightcurve due to variations on time-scales of 
to ly^^ (e.g. van der Klis 1997). Thus the total fractional 
RMS variability of the lightcurve is given by the square 
root of the integral of the power spectrum across all mea- 
sured frequencies, Vmin to i^Nyq- Under this normalisation, 
the constant level of power contributed to all frequencies by 
the Poisson noise in the lightcurve is equal to 2{fi + B) / fi^ , 
where B is the mean background count rate. Using this nor- 
malisation allows us to compare power spectra measured by 
different instruments and power spectra of different sources, 
and take account of the linear RMS-flux relation recently 
discovered in AGN and XRBs (Uttley & M'^Hardy 2001, 
and see Section 3.4). 

For each source in our sample we have lightcurves for three 
observing schemes, which we use to measure power spectra 
over three different frequency ranges to produce the broad- 
band power spectrum: 

1. A long-term monitoring lightcurve incorporating all 
monitoring data, to measure the low-frequency power spec- 
trum (~ 10"* Hz-10"'^ Hz). 

2. An intensive monitoring lightcurve, to measure a 
medium-frequency power spectrum (~ 10"'' Hz-10"^ Hz). 

3. A long-look lightcurve (two such lightcurves for 
NGC 3516) to measure the high-frequency power spectrum 
~ lO"'^ Hz-lO"* Hz. 

Additionally, for the most variable sources MCG-6-30-15 
and NGC 5506, which show significant variability on time- 
scales less than 1 ks, we measure a very-high-frequency 
(VHF) power spectrum (~ 4 x 10~* Hz-10~^ Hz) using con- 
tinuous ~ 2.5 ks segments of the PCA lightcurves (i.e. be- 
tween Earth-occultations of the source), binned to 16 s reso- 
lution. We do not include VHF power spectra for NGC 5548 
and NGC 3516, since they show no significant source power. 



other than the small amount expected at the lowest frequen- 
cies due to red-noise leakage of variations which are sam pled 
by the high-frequency power spectrum (see Section 3.3). 
In order to minimise any distortion, the power spectra are 
made from lightcurves binned up to the maximum sam- 
pling interval of the observing scheme under consideration 
(i.e. the lightcurve resolution is two weeks or 1209.6 ks for 
the long-term monitoring lightcurves, 86.4 ks for the inten- 
sive monitoring lightcurves, except for NGC 3516 where we 
bin the long-term and intensive monitoring lightcurves to 
4.3 days and 12.8 hours respectively). Long-look lightcurves 
•are binned to 2048 s for the purposes of making the high- 
frequency power spectra, so that gaps due to Earth occulta- 
tion are minimised to be no more than one bin wide. Empty 
lightcurve bins in the binned-up monitoring and long-look 
lightcurves are filled by linearly interpolating between adja- 
cent filled bins. No rebinning or interpolation was applied to 
the 16 s lightcurves used to determine the VHF power spec- 
tra of NGC 5506 and MCG-6-30-15, since only continuous 
sections of the lightcurves were used to estimate the power 
spectrum. 

The total lightcurve durations, bin widths, mean fiuxes, frac- 
tional RMS variability (after subtracting the Poisson noise 
contribution to variance) and power-spectral Poisson noise 
levels for each lightcurve are given in Table |^. Note that 
mean flux and fractional RMS are calculated based on the 
quoted bin widths, i.e. the contributions to mean fiux and 
fractional RMS from each bin are equally weighted, so that 
bins containing many data points (e.g. 2-week wide bins 
which contain daily or twice-daily observations) do not con- 
tribute more to the mean flux or variance than bins which 
contain a single data point. 

We measured each power spectrum using the method and 
normalisation outlined above. In order to reduce the scat- 
ter in power-spectral points, which fiuctuate wildly for a 
stochastic process such as red or white-noise, we binned 
the logarithm of power at each frequency (see Papadakis 
& Lawrence 1993) in logarithmically spaced frequency bins, 
separated by a factor of 1.3 in frequency but with a minimum 
of two measured powers per bin, so that the bin spacing is 
larger at the lowest frequencies sampled by each power spec- 
trum. The VHF power spectra for NGC 5506 and MCG-6- 
30-15 were calculated by measuring separate power spectra 
for each continuous lightcurve segment, averaging them and 
binning in logarithmically spaced bins separated by a factor 
of 1.3 in frequency. The resulting broadband power spectum 
for each object is shown in Figure ^ 

Inspection of the power spectra in Figure ^ shows that they 
do flatten at low frequencies. However, we cannot immedi- 
ately assume that this flattening is real and representative of 
the shape of the true, 'underlying' power spectrum, for the 
following reasons: 1) First of all, we cannot estimate reliable 
errors for all but the VHF power spectra, especially for the 
points at the low-frequency end of each power spectrum, due 
to the small number of points which contribute to each fre- 
quency bin, 2) although rebinning and interpolation result in 
evenly sampled lightcurves, they also introduce distortions 
in the estimated power spectra which are difficult to predict 
a priori. Furthermore, even if these distortions are mini- 
mal, the estimation of red-noise power spectra is affected by 
potentially serious distortions due to aliasing and red-noise 
leak, which are dependent on the original sampling pattern. 
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Table 2. Parameters of lightcurves for use in PSRESP 
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T and ATsamp are the lightcurve duration and sampling interval (in seconds), ^ and o-frac are the lightcurve mean flux (in count s~^) 
and fractional RMS respectively and Pnoise is the Poisson noise level expected in the power spectrum due to counting statistics (in 
fractional RMS-squared units, Hz~^). Notes: " Details given in the table are for the iJJfOSylT lightcurve, the RXTE lightcurve used to 
measure the power spectrum at the highest frequencies has /i = 28.1, ^*noisc — 0.10. Lightcurve obtained 22-26 May 1997. 
Lightcurve obtained 13-16 April 1998. 



Finally, 3) we must consider the possibility that the underly- 
ing power spectra are not stationary, but vary on time-scales 
comparable to the length of our campaign, so that it is not 
valid to combine power spectra taken at diflterent times and 
over different intervals. We consider these problems in more 
detail in the remainder of this section. 



3.2 Error estimation 

The smooth functions used to fit the power spectra of noise 
processes such as red-noise, white-noise and the composite 
broadband noise represent the average power spectrum of 
the underlying noise process, Pproc(i')- However the hght 
curve which we measure is a stochastic realisation of that 
process and results in an observed power spectrum Pobs(i') 
which fluctuates randomly about Pproc(i^), following a 
distribution with 2 degrees of freedom and standard devia- 
tion at any frequency v equal to Pprociy) (e.g. Timmer & 
Konig 1995). Therefore, in order to recover the underlying 
power spectrum of the process Parody) directly from the 
data, we must determine the mean power spectrum by av- 
eraging many observed power spectra, using the spread in 
the observed power measured at each frequency to estimate 
the standard error on the mean. Because the distribu- 
tion is exponential the power at a given frequency fluctu- 
ates wildly, so the number of power spectra averaged must 
be large (> 50), in order that the standard error is reliable. 
This problem has been discussed extensively by Papadakis 
& Lawrence (1993), who show how more reliable estimates 
of smoothed power (and the standard error) can be obtained 
by averaging fewer power spectra (~ 20) if we instead aver- 
age the logarithm of power rather than the power. 
Unfortunately, we cannot use this method of error estima- 
tion to constrain the shape of the power spectra we measure 
here (other than the VHF power spectra) , because there are 
not enough data points to estimate reliable errors, especially 



at the lowest frequencies measured in each power spectrum. 
Therefore we must discard the requirement that the data 
are used to directly estimate a reliable mean power-spectral 
shape and errors and instead use a Monte Carlo technique, 
using simulated lightcurves to estimate the power-spectral 
shape and uncertainty for a range of specified models and 
test these against the data. Using this approach we can esti- 
mate reliable uncertainties even in the limit of small-number 
statistics and hence use the full range of power-spectral fre- 
quencies available to us. Furthermore we can take account of 
the distorting effects of lightcurve sampling, which we detail 
below. 



3.3 Power-spectral distortion due to lightcurve 
sampling 



As we described in Section 3.1 , the lightcurves that we use to 
calculate the power spectra are rebinned to an even pattern 
and any empty bins are filled with interpolated flux measure- 
ments. The use of the Monte Carlo technique mentioned in 
the previous section can take account of the distorting effects 
of rebinning and interpolation on the power spectrum (see 



Section 4.1). However, the estimation of a red-noise power 
spectrum, even from an evenly sampled lightcurve, is not 
free from distortions. 

Consider an underlying continuous lightcurve r{t) whose 
Fourier transform is R{u), on which we impose a sampling 
pattern w(t) so that io(t) = 1 when we sample r{t) and zero 
otherwise. The resulting observed lightcurve, f{t) is given 
by: 

f{t) = rit)w{t). 

Applying the Convolution theorem of Fourier transforms, 
the Fourier transform of /(t), F(y) is then given by the con- 
volution of R{v) and the Fourier transform of the sampling 
pattern Wiy) (known as the 'window function'), i.e. 
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Frequency (Hz) 

Figure 5. Raw broadband power spectra of MCG-6-30-15, NGC 5506, NGC 5548 and NGC 3516. The dashed line shows the low-frequency 
part of the power spectrum, made from the total monitoring lightcurves, while the dotted line shows the high-frequency part made from 
long-look lightcurves. Solid lines mark the medium- frequency power spectrum (made from the intensive monitoring lightcurves), VHF 
power spectra for MCG-6-30-15 and NGC 5506, and the power spectrum of the second long-look observation for NGC 3516. Poisson 
noise levels have not been subtracted from the power spectra. 



F{u) = R{u) *W{u). 

Therefore the Fourier transform (and hence the power spec- 
trum) of the observed lightcurve is distorted from the true 
underlying power spectrum by the sampling pattern im- 
posed on the underlying lightcurve. Qualitatively we can 
distinguish two significant components to this distortion in 
the case of red-noise power spectra, red-noise leak and alias- 
ing. 

Significant power below the minimum frequency sampled 
by the power-spectrum {vmin = T^^) causes long-time-scale 
trends in the lightcurve which cannot be distinguished from 
smaller amplitude trends on the time-scales which are sam- 
pled by the power spectrum. The result of this red-noise leak 
is that additional power is transferred across the entire mea- 
sured power spectrum, with an amplitude dependent on the 
amount of power at frequencies below !/min (and hence the 
amount of red-noise leak is model dependent and stochastic). 
Fortunately, the effects of red-noise leak can be accounted for 
using the Monte Carlo technique mentioned earlier, by en- 



suring that the simulated 'underlying' lightcurves are much 
longer than the observed lightcurves (see Section |4t| ). 
If a lightcurve is not continuously sampled (i.e. it is sampled 
for a duration ATbin at sampling intervals ATsamp, where 
ATbin ^ ATsamp) then variations on time-scales shorter 
than ATsamp (i.e. corresponding to power above the Nyquist 
frequency, !^Nyq) cannot be distinguished from (and therefore 
appear to contribute to) variations on longer time-scales. 
The result is that power is shifted or 'aliased' to lower fre- 
quencies from frequencies above v^yq- Technically, the ef- 
fect of aliasing is to transfer the power at a frequency above 
the Nyquist frequency, ^'Nyq + Ai^ to a frequency below the 
Nyquist frequency h^yq—Ai/, i.e. the power is refiected about 
I'Nyq (e.g. see van der Klis 1989). Hence, the amount and 
form of aliasing in the observed power spectrum is depen- 
dent on the underlying power-spectral shape and amplitude. 
However, since for all but the steepest broadband-noise type 
power spectra, the power at 2!/Nyq is not much less than the 
power at i/Nyq, the result of aliasing can be approximated 
by adding a constant level of power to the underlying power 
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spectrum. For lightcurves with initial time resolution ATbin 
(prior to any rebinning, e.g. 1 ks in the case of our moni- 
toring lightcurves), we expect that variations with frequen- 
cies higher than ~ (2ATbin)~^ will be smoothed out and 
will not contribute significantly to aliasing. For this reason, 
we expect the total integrated power transferred to the ob- 
served power spectrum by aliasing to be roughly equal to 
the integrated power between the Nyquist frequency and 
(2ATbin)~^. As a first approximation, we assume that this 
power will be distributed evenly to all sampled frequencies, 
with the constant power, Pc, added to all frequencies be- 
cause of aliasing given by: 

^(2ATbi„)-l 

Pc = / P{p) diy. 

l^Nyq - i^min J ^^^^ 

The effects of aliasing on power spectra of different under- 
lying slopes is shown in Figure ^. The average aliased power 
spectra were constructed from fOOO simulated lightcurves, 
each lightcurve corresponding to 2 years of evenly spaced 
weekly 1 ks snapshots. The assumed underlying power spec- 
trum (shown by the dashed lines) was cut off above 10~* Hz, 
to reduce the red-noise leak contribution to the power spec- 
tra. We also plot our estimate of the aliased power spec- 
trum obtained using the constant-power approximation out- 
lined above (dotted lines). Note that the flattening of the 
power-spectral slope is more pronounced for flatter under- 
lying slopes, as expected. At high frequencies, the agree- 
ment between the aliased power spectrum and our estimate 
is very good, except for the steepest power spectrum where 
the amount of aliasing is small anyway. 
Fortunately, the distorting effects of aliasing can be properly 
taken account of using the Monte Carlo technique mentioned 
earlier. The approximation to the distorting effect of aliasing 
presented here will prove useful for reducing the resolution 
of the simulated lightcurves required by the technique, as 
we will show in Section 4.1, 



3.4 Stationarity of the power spectrum 

A key assumption we must make, in order to combine 
lightcurves obtained at different times to measure the broad- 
band power spectrum, is that the underlying power spec- 
trum is stationary, i.e. its amplitude and shape do not 
change over time (note that power-spectral non-stationarity 
is not the same as lightcurve non-stationarity, which is ex- 
pected on short time-scales for red-noise processes). If the 
underlying power spectrum were non-stationary, so that it 
happened to be steeper while the high-frequency power spec- 
trum was measured than when the medium or low frequency 
power spectra were measured, then these changes in power- 
spectral slope would masquerade as a flattening in the power 
spectrum. 

The power spectra of BHXRBs are known to be non- 
stationary on a variety of time-scales, showing drastic 
changes in shape between the well-known low and high states 
(e.g. Cui et al. 1997). Transitions between these states occur 
on time-scales of months to years and less drastic changes 
in power-spectral shape within the low and high states oc- 
cur on time-scales of days or weeks (e.g. changes in the 
high-frequency slope and the position of the power-spectral 
breaks, Belloni & Hasinger 1990, Cui et al. 1997). If the 
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Figure 6. The effects of aliasing on power-spectral shape. The 
grey dashed line represents the true underlying power spectrum, 
while the solid black line shows the average observed power spec- 
trum after distortion by the effects of discrete sampling at weekly 
intervals, obtained from averaging the power spectra of 1000 sim- 
ulated lightcurves of 2- year duration (see text). The dotted line 
shows our approximation to the average observed power spec- 
trum, obtained by assuming that the effect of aliasing is to add 
a constant to the underlying power spectrum. Underlying power 
spectral slopes from top to bottom are a =2.0, 1.5 and 1.0 re- 
spectively. 



variability time-scales for this kind of non-stationarity scale 
linearly with black hole mass, we would expect to see sim- 
ilar changes in the power spectra of AGN on time-scales of 
centuries or longer - much greater than the time-scale we 
can sample. This picture is supported by the fact that, to 
date, no hard evidence for non-stationarity in power-spectral 
shape has been reported for AGN. Hence it is probably safe 
to assume that the shapes of our target's power spectra do 
not vary during the course of our monitoring campaign. 
An alternative concern is that the amplitude of the power 
spectrum varies over time, even though the shape does not. 
In a separate paper, we report the intriguing (and unex- 
pected) result that the X-ray RMS variability of Cyg X- 
1 and the accreting millisecond pulsar SAX J1808. 4-3658 
scale linearly with the local mean flux (Uttley & M'^Hardy 
2001). In other words, the fractional RMS (i.e. RMS di- 
vided by mean flux) measured for any fixed-length segment 
of the lightcurve is constant (subject to the stochastic vari- 
ability inherent in the lightcurves), irrespective of the seg- 
ment's mean flux. We also confirm that the same property 
applies to the lightcurves of AGN and hence appears to be 
an important characteristic of the broadband noise variabil- 
ity in compact accreting systems. It is obvious then that 
we must always normalise our lightcurves by their mean be- 
fore calculating their power spectra or, as we do here, divide 
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their power spectra by the square of the mean, otherwise 
the shape of our broadband power spectrum will be depen- 
dent on the mean object flux at the time we measured each 
lightcurve. 



4 PSRESP: A RELIABLE METHOD OF 

ESTIMATING POWER-SPECTRAL SHAPE 

In order to take proper account of the errors on the power 
spectrum and the distorting effects due to rebinning, red- 
noise leak and aliasing, we must use a Monte Carlo tech- 
nique to estimate the underlying power-spectral shape (and 
its uncertainty). The basic concept is to simulate a large 
number of continuous lightcurves with a known underlying 
power-spectral shape, apply the sampling pattern of the ob- 
served lightcurve, rebin and interpolate as necessary, ob- 
tain the resulting distorted power spectra and average them 
to determine the mean shape of the distorted model power 
spectrum. The spread of the simulated power spectra about 
the mean can be used to estimate the errors on the observed 
power spectrum. Applying the sampling pattern to a simu- 
lated (i.e. model) lightcurve before measuring its power spec- 
trum is equivalent to convolving the window function of the 
sampling pattern with the Fourier transform of the model 
lightcurve to obtain the Fourier transform of the observed 
lightcurve. The effect is then similar to that used in X-ray 
spectral measurement, when a model spectrum is postulated 
and then convolved with the 'response function' of the de- 
tector to yield the 'observed' model spectrum which can be 
tested against the data. Hence, this Monte-Carlo technique 
for estimating the true power-spectral shape is known as the 
'response method'. 

The response method was applied by Done et al. (1992), to 
measuring the power spectrum of unevenly sampled Ginga 
lightcurves and has subsequently been applied to ROSAT 
data (Green, M'^Hardy & Done 1999). In both cases the 
method was applied only to a single lightcurve, of relatively 
short (a few days) duration. Here, however, we wish to fit the 
same power-spectral model to several power spectra in or- 
der to estimate the shape of the broadband power spectrum. 
We should not fit power spectra measured over different fre- 
quency ranges separately, because the distorting effects of 
aliasing and red-noise leak, which we must take account of, 
are dependent on the shape of the underlying power spec- 
tral shape outside the frequency range measured by a single 
power spectrum. Furthermore, the response method of Done 
et al. (1992) estimates a goodness of fit based on the stan- 
dard statistic, which is not a reliable estimator for the 
distorted power spectra considered here and hence cannot be 
used to formally estimate fit probabilities. Therefore, for our 
purpose of reliably constraining the shape of the broadband 
power spectrum of AGN, we have developed a more sophis- 
ticated technique based on the response method, which we 
call PSRESP (for Power Spectral RESPonse). In the remainder 
of this section, we describe PSRESP in detail. 



4.1 Lightcurve simulation 

We simulate lightcurves using the method of Timmer & 
Konig (1995), which is superior to the commonly used 
method of summing sine waves with randomised phases in 



that the power-spectral amplitude at each frequency is ran- 
domly drawn from a distribution, as should be the case 
for noise processes, and not fixed at the amplitude of the 
underlying power spectrum. We create the fake 'observed' 
lightcurves as follows. 

First, we specify a power-spectral model which we wish to 
test against the data (which is some continuous function 
such as a power-law, with or without breaks) . The normali- 
sation of the model power spectrum is a multiplicative factor 
which is carried through any convolution with the window 
function (i.e. only the power-spectral shape is distorted by 
sampling) . Hence we can choose an arbitrary model normal- 
isation and simply renormalise the resulting distorted model 
power spectrum to obtain the best possible fit to the data. 
For each observed lightcurve used to measure the broadband 
power spectrum, we simulate A*" continuous lightcurves using 
the given power-spectral model (where A'' is large, between 
100 and 1000). Ideally, the time resolution of the simulated 
lightcurves, ATsim should be the same as the initial resolu- 
tion of the observed lightcurve ATbin, which for our mon- 
itoring data is the typical exposure time of the snapshot 
observations. This is because any variations on time-scales 
down to this resolution will contribute to aliasing. Although 
we have derived an analytical approximation to the effect of 
aliasing on the power spectrum, this only tells us the average 
effect of aliasing for lightcurves which are evenly sampled. 
Due to the stochastic nature of the lightcurves, the actual 
power above the Nyquist frequency is variable and hence 
aliasing contributes a variable amount to the power spec- 
trum which adds to the uncertainty in its shape. However, 
for the monitoring lightcurves, which consist of only 1 ksec 
snapshots but may be months or years long, the requirement 
that the simulated lightcurves also have 1 ksec resolution in- 
creases the computation time prohibitively. 
In practice, we can limit the resolution of the simulated 
lightcurves to be ATaim < O.lATsamp. This is because, 
for lightcurves whose power spectra are steep at high fre- 
quencies, like those we measure here, the uncertainty in 
the amount of aliasing is dominated by the uncertainty in 
the large amount of power at frequencies not much greater 
than the Nyquist frequency. We can then estimate the much 
smaller contribution to aliasing due to power at frequen- 
cies greater than sampled by the simulated lightcurves using 
our analytical approximation, i.e. calculating the integrated 
power of the model power spec trum between (2Arsini)~^ 
and (2Arbin)~^ (see Section 3.3) and incorporating the re- 
sulting constant values into the final simulated power spec- 
tra. 

In order to take account of red-noise leak, the simulated 
lightcurves must be significantly longer (e.g. by a factor 10 
or more) than the observed lightcurve, so that there is power 
at frequencies lower than the minimum frequency sampled 
by the observed lightcurve. We can minimise the cost of 
simulating excessively long lightcurves by making a single, 
very long lightcurve of length NT, which is subsequently 
split into the A'' desired segments. Note that although the 
longest time-scale trends in the total simulated lightcurve 
contribute the same amount to the red-noise leak in the 
power spectra of all simulated lightcurve segments, the 
bulk of red-noise leak (and the power-spectral uncertainty 
it introduces) is due to variations on shorter time-scales close 
to the observed lightcurve duration and hence will be sta- 
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tistically independent between segments. 
Once a continuous lightcurve is simulated, it is resam- 
pled using the same sampling pattern as the observed 
lightcurve. The resampled lightcurve is then rebinned and 
empty bins interpolated in the same manner as for the 
observed lightcurve. The power spectrum of the resulting 
lightcurve is then measured. 



4.2 Determining the goodness of fit of the model 

The power spectrum of each simulated lightcurve is calcu- 
lated and binned in the same way as the power spectrum of 
the original observed lightcurve. The model average power 
spectrum (corresponding to the given model and sampling 
pattern of the observed lightcurve) is then given by the mean 
of the A'^ simulated power spectra, Psim(i'). For each fre- 
quency v sampled by the binned power spectrum, we deter- 
mine the RMS spread of simulated model powers about the 
mean and define this spread to be the 'RMS error' on the 
power at that frequency, APsini(i^). We next define a statis- 
tic, which we call Xdisti which is calculated from the model 
average (and RMS error) and observed power spectra of each 
lightcurve: 



Xdist 



where Umin and //max are respectively the minimum and max- 
imum frequencies measured by Pobs(i'). We measure xlist 
for each input power spectrum (i.e. low, medium and high- 
frequency), and sum to yield a total Xdist fo'' the model 
with respect to the data. Note that although the approach 
of assigning error bars to the model rather than the data is 
unusual, it is technically valid and strictly speaking the more 
correct thing to do, since the statistic is defined on the 
basis of the variance of the model population, which error 
bars conventionally determined from the data are supposed 
to approximate (e.g. see discussion in Done et al. 1992, Sec- 
tion 4.2). 

Next, we find the best-fitting normalisation of the power- 
spectral model by renormalising the Psim(i/) for each input 
power spectrum by the same factor k, varying k until the 
total Xdist is minimised. 

Having determined the minimum Xdist ioi the model com- 
pared to the data, we next estimate what goodness of fit this 
value of Xdist corresponds to. The Xdist is not the same as the 
standard x^ distribution, because the Pobs(i')s are not nor- 
mal variables (since the number of power spectrum estimates 
averaged in each frequency bin is small) . Therefore we must 
estimate a reliable goodness of fit using the distribution of 
Xdist s-nd not the well-known x^ distribution. For each in- 
put low, medium and high-frequency power spectra, we have 
already created corresponding simulated power spectra 
which can be used to determine the distribution of Xdist for 
that particular model and lightcurve sampling pattern. In 
order to determine the distribution of Xdist j we should cal- 
culate the values of Xdist for all possible combinations of 
simulated low, medium and high-frequency power spectra. 
However, since the number of such combinations (A'^^) may 
be extremely large, we reduce the number of combinations 
we sample to M (where M > 1000), randomly selected to 
ensure an accurate estimate of the distribution of Xdist- 



We sort the M simulated measurements of the total Xdist 
into ascending order. The probability that the model can be 
rejected is then given by the percentile of the simulated xlist 
distribution above which Xdist exceeds that measured for the 
observed power spectra. Note that this method of using sim- 
ulated data sets to estimate goodness of fit in the absence of 
a well-understood fit statistic is well-known and described 
in Press et al. (1992), Section 15.6. 

4.3 Incorporating VHF power spectra 



As described in Section 3.1, we use continuous segments of 
the 16 s binned long-look lightcurves to make VHF power 
spectra for our most variable targets, by averaging the power 
spectra of all the segments and determining the standard 
error from the spread in power at each frequency. The stan- 
dard errors estimated from the data are reliable, since we 
average > 20 power spectra and bin the logarithm of power 
(according to the method of Papadakis & Lawrence 1993 and 
see Section 3.S). We can therefore use the measured VHF 



power spectra and their errors as they are, without hav- 
ing to estimate errors using simulations of high-resolution 
lightcurves which would be very computationally intensive. 
However, if we simply compare the VHF power spectrum 
with the underlying undistorted model shape, we ignore the 
effects of red-noise leak which could be significant in dis- 
torting the shape of the VHF power spectrum, especially if 
the underlying power spectrum docs not flatten significantly 
until far below the minimum frequency sampled by the VHF 
power spectrum. Therefore, we need to take account of the 
effects of red-noise leak on the model power-spectral shape 
at high frequencies. Note that, because the segments are con- 
tinuous and binned into high-resolution time bins of width 
16 s, the effects of aliasing are not important in this case. 
To determine the effects of red-noise leak on the VHF power 
spectrum, we simulate 1000 lightcurves, each made to be at 
least 64 times longer than ly^l^ (where i/min is the minimum 
frequency sampled by the VHF power spectrum), with res- 
olution ATsim smaller than |i^max where i/max is the maxi- 
mum frequency which contains significant power above the 
noise level (typically around 3 x 10~^ Hz) and is chosen so 
that the ratio of i^max to Vmin is a power of 2. Power spectra 
of the lightcurves sampled to have duration v^l^ may then 
be made using the Fast Fourier Transform, which allows 
a VHF model average power spectrum for 1000 simulated 
lightcurves to be determined very rapidly. The VHF model 
average power spectrum is then used in place of the underly- 
ing model power spectrum, while the errors determined from 
the observed VHF power are used as errors on the model. 
The contribution of the VHF power spectrum to the total 
Xdist is determined by comparing the observed power spec- 
trum with the simulated model average power spectrum, us- 
ing the standard errors estimated from the data. The con- 
tribution of the VHF power to the goodness of fit of the 
model is obtained as follows: for each of the M random com- 
binations of simulated power spectra used to estimate the 
goodness of fit, a VHF power spectrum is simulated by ran- 
domly selecting the power at each measured frequency, from 
a Gaussian distribution of mean equal to the model average 
power and standard deviation equal to the standard error at 
that frequency. The Xdist of the simulated VHF power spec- 
trum is determined and included in the total Xdist measured 
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for that selection of simulated power spectra. The goodness 
of fit of the model is then estimated as described in the 
preceding section. 



4.4 Summary of the psresp method 
We summarise the PSRESP method as follows: 

1. Measure the observed power spectrum of each (re- 
binned and mean-subtracted) observed lightcur ve a nd bin 
up the power spectrum as desired (see Section 3.1). Mea- 
sure the VHF power spectrum if require d, d etermining its 
errors directly from the data (see Section 4.3). 

2. Specify the underlying power spectral model to be 
tested against the data. For the given set of parameters, 
simulate A'' lightcurves which are realisations of the model 
and apply the sampling pattern of the o bse rved lightcurve 
to the simulated lightcurves (see Section [4.l[ ). 

3. Calculate the power spectrum of each re-sampled sim- 
ulated lightcurve using the same method used to measure 
the observed power spectrum. Determine the model aver- 
age power spectrum from the simulated power spectra, 
and assign error bars equal to the RMS spread in simulated 
power at each frequency (see Section 4.2). 



The two steps above should be repeated for each lightcurve 
(i.e. long-term, intensive and long-look), to make simulated 
model average power spectra corresponding to each observed 
power spectrum. The model average VHF power spectrum 
should also be determined at this point (if required) , accord- 



ing to the method outlined in Section 4.3 

4. Estimate the statistic xlist (summed over all input 
power spectra) for the observed versus the model average 
power spectrum, and vary the normalisation of the model 
to minimise Xdist obtain the best-fitting normalisation 
(include the VHF power spectrum in this determination if 
required, using standard errors determined from the data, 
see Section 4.3). 

5. Determine the xlist for M randomly selected combi- 
nations of the simulated power spectra. If a VHF power 
spectrum is included, measure its contribution to each sim- 
ulated total Xdist from a random realisation of the model 
average VHF powe r sp ectrum, according to the method de- 



scribed in Section 4.3. Sort the simulated distribution of 
Xdist into increasing numerical order - the model is rejected 
at a confidence equal to the percentile of the simulated xlist 
distribution above which Xdist exceeds that measured for the 
observed power spectra (see Section ^!^) . 

Using the method described above, any given model can 
be tested against the data. By stepping through a range of 
model parameters and repeating steps 2-5, large regions of 
the model parameter space may be searched and confidence 
regions may be determined. 



5 RESULTS 

We will now apply the PSRESP method described in Section ^ 
to the lightcurves of our sample, in order to determine if the 
broadband power spectra of Seyfert galaxies really flatten 
towards low frequencies and to try to constrain models for 
any flattening which we see. 



5.1 Do the broadband power spectra really 
flatten? 

To determine if the power spectra flatten, we will test a 
simple power-law model for the underlying power spectrum, 
fmod(i^) of the form: 



-Pmod(l^) = A 



where A is the amplitude of the model power spectrum at 
a frequency vo, a is the power-spectral slope and Cnoiso is a 
constant value which is flxed at the Poisson noise level for 
the lightcurve. Note that the Poisson noise level is included 
in the model rather than subtracted from the power spectra 
before model fitting, because the power spectra are binned 
logarithmically (so constants in linear space may not simply 
be subtracted). This is particularly important for the VHF 
power spectra of NGC 5506 and MCG-6-30-15, whose stan- 
dard errors are determined in logarithmic space, and also for 
high-frequency power spectra in general, which are close to 
the Poisson noise level, since fluctuations in the power spec- 
trum lead to some measured powers lying below the Poisson 
noise level (so subtraction of this level would lead to nega- 
tive measured powers). 

The model is tested against the measured power spectra 
by stepping through values of a from 1.0 to 2.4 in incre- 
ments of 0.1 (i.e. test the model with a = 1.0, 1.1, 1.2 etc.). 
These values of a cover the range of reasonable values which 
could possibly be fltted to the data. Probabilities that the 
measured power spectra could be a realisation of the model 
are calculated by PSRESP, as described in Section 4.2 using 
A*' — 1000 simulated lightcurves to determine the distorted 
model average power spectrum. The distribution of Xdist of 
the realisations of the model is determined for each value of 
a by measuring xlist for ^ = 10* randomly selected sets 
of simulated power spectra. The simulated lightcurves have 
time resolutions ATsim given in Table ^. Additional distor- 
tion in the power spectrum due to model power at frequen- 
cies greater than (2Arsini)~^ is calc ulated directly from the 



model, as described in Section 4.1. Distorted VHF model 



power spectra, which take account of red-noise leak in the 
VHF power spectra included in the broadband power spec- 
tra of MCG-6-30-15 and NGC 5 506 , are determined using 
the method described in Section 4.3. 



The best-fitting values of a, and corresponding confldences 
that the single power-law model is rejected by the data, are 
given in Table H. The flrst and second of each of these values 
shown for NGC 5506 correspond to fits without or includ- 
ing the EXOSAT data respectively. The simple power- law 
model, without any flattening is rejected at better than 99% 
confidence for MCG-6-30-15, better than 90% confidence (or 
close to 99% confidence including the EXOSAT data) for 
NGC 5506 and better than 95% confidence for NGC 3516. 
Only for NGC 5548 is the model not rejected at a significant 
confidence. The best-fitting models are compared with the 
measured power spectra in Figure ^ 

It is apparent from these plots that the simple power-law 
model does not fit the observed power spectra of MCG-6- 
30-15, NGC 5506 and NGC 3516, even after allowing for the 
distorting effects of sampling, because the intrinsic power 
spectrum of each of these objects does indeed fiatten towards 
low frequencies. There is no significant evidence for flatten- 
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Table 3. Time resolution of simulated lightcurves used in PSRESP. 



MCG-6-30-15 
NGC 5506 
NGC 5548 
NGC 3516 



Long-term 
86400 s 
86400 s 
86400 s 
46080 s 



Intensive 
8640 s 
8640 s 
8640 s 
4608 s 



Long-look 
512 s 
512 s 
512 s 
512 s 



10' 
lo" 
10^ 
10" 
1000 
100 
10 

1 
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Figure 7. Comparison of best-fitting model average power spectra with the observed power spectra for the single power-law model 
described in the text. Open squares mark the low-frequency model average, simple crosses mark the high-frequency model average while 
filled squares mark the medium-frequency model average, the VHF model average for MCG-6-30-15 and NGC 5506, and the high- 
frequency model average for the second long-look observation of NGC 351 6. Note that the error bars represent the RMS error in the 
simulated power spectra used to calculate Xjist described in Section [4.2| . 



ing at low frequencies in the power spectrum of NGC 5548. 
Figure ^ shows the fit probability measured at each input 
value of Q for NGC 5548, which demonstrates how PSRESP 
is capable of finding well-defined probability maxima in the 
same way that fitting can, using more conventional data 
sets. 

The simple power law model with no frequency breaks can 
be rejected at better than 95% confidence for all objects ex- 
cept NGC 5548. The next step is to try to fit the observed 
power spectra with more complex models which fiatten at 
low frequencies, in particular, can we distinguish between 
models where the power spectrum fiattens to a = or a = 1, 



Table 4. Results from fitting broadband power spectra of four 
Seyfert galaxies with a simple unbroken power-law model. 



MCG-6-30-15 
NGC 5506 
NGC 5548 
NGC 3516 



Best-fitting a 
1.5 
1.4/1.5 
1.6 
1.8 



Rejection confidence 
99.8% 
90.6%/98.6% 
67% 
96.6% 



Quoted values for NGC 5506 correspond to fits 
excluding/including the EXOSAT data,. 

and can we constrain any characteristic frequencies for the 
fiattening? 
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Figure 8. Probability that a single power-law of slope a is ac- 
ceptable to describe the broadband power spectrum of NGC 5548. 



Figure 9. The form of the knee model power spectrum. 



5.2 Fitting simple models for the power-spectral 
flattening 

Although we can say with confidence that the power spec- 
tra of three of the objects in our sample flatten towards low 
frequencies, it is not clear what form this flattening takes. 
In this work, we will restrict ourselves to considering two 
simple models, a 'knee' model based on the low-frequency 
flattening seen in the power spectra of BHXRBs in the low 
state (i.e. the low-frequency break in Cyg X-1) and a 'high- 
frequency break' model which assumes that the flattening is 
due to a frequency break in the power spectrum analogous 
to the high-frequency break seen in the power-spectrum of 
Cyg X-1. Under the knee model, the power-spectrum turns 
over to a slope a — below some 'knee frequency', whereas 
under the high-frequency break model, the power spectrum 
breaks sharply to a = 1 below some 'break frequency'. More 
complex models, consisting of multiple frequency-breaks and 
a variety of power-spectral slopes, or a number of broad 
Lorentzians (e.g. Nowak 2000) might also successfully ex- 
plain the data. However, computational constraints limit the 
testing of a large variety of models for the flattening we see 
and moreover, as we shall discover, the data do not yet war- 
rant these kinds of models as the observed power spectra 
can be fltted adequately by the simple models we test here. 
We will fit these two simple models for the fiattening to the 
power spectra of all the Seyfert galaxies in our sample, in- 
cluding NGC 5548 so as to set upper limits on any knee or 
break frequencies. 



5.2.1 The knee model 

We first test the knee model for the power spectrum, which 
has the form 



A 



1 + 



j/2 ' 



where A is the constant amplitude of the power-spectrum 
at zero slope, i^knoc is the 'knee frequency' and a is now de- 
fined as the power-spectral slope above the knee frequency. 



The shape which this function describes can be seen in Fig- 
ure ^ We now test this model against the measured broad- 
band power spectra, to see if it can explain the fiattening 
we see. Using the equation given above for the underlying 
power spectral shape (also including the constant Poisson 
noise level), we can fit the model in the same way as fit- 
ting a simple power law in the previous section. The free 
parameters to be stepped through are a, which is again in- 
cremented in steps of 0.1 between 1.0 and 2.4, and fknoc 
which is stepped through by multiplicative factors of 2, from 
10"* Hz to 10"^ Hz, since a very broad range in frequency 
must be covered. Approximately 200 pairs of a and J^knoo 
must be tested (as opposed to only 15 parameters when fit- 
ting the simple power law in the previous section), so to 
save on computing time, the number of lightcurve simula- 
tions used to estimate each model average power spectrum 
for each pair of parameters is reduced from A'^ — 1000 to 
N = 100. The xi ist distribution is obtained by determining 
Xdist foi' ^ ~ 1000 sets of simulated power spectra. 
The best-fitting parameters and probabilities are shown in 
Table ^. Contour plots for each of the knee model fits (not 
including the fit of the NGC 5506 broadband power spec- 
trum which excludes the EXOSAT data), together with the 
best-fitting model average power spectra, are shown in Fig- 
ure |l^. The contour plots show that the acceptable regions 
are broad and well-defined, rather than consisting of very 
many separate 'islands', which implies that using only 100 
simulated lightcurves per measured power spectrum is suffi- 
cient to determine reliable maxima in the probability space. 
As Table ^shows, the knee model fits the power spectra of all 
the objects in our sample adequately. In all cases, the model 
can be accepted at a confidence level better than 10%. Note 
that the broadband power spectrum of NGC 5548 is consis- 
tent with this model, although the lower confidence limits 
on the knee frequency cannot be defined, in agreement with 
the acceptable simple power law fit to these data. The fact 
that the knee model provides a good fit to the broadband 
power spectrum of NGC 5506 including the EXOSAT ded,a,, 
is consistent with the power spectrum of NGC 5506 being 
stationary over time-scales as long as a decade. 
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Table 5. Results from fitting broadband power spectra of four Seyfert galaxies with a knee model 



MCG-6-30-15 
NGC 5506" 
NGC 5506'' 
NGC 5548 
NGC 3516 



Rejection confidence 
81% 
42% 
32% 
53% 
83% 



!^knco/10-«' Hz 

5.12 (2.56-10.24) 
0.64 (0.0-10.24) 
2.56 (0.16-10.24) 
0.02 (0.0-1.28) 
2.1 ±0.3 0.64 (0.32-1.28) 



a 

1.8 ±0.1 

1 T+O'^ 
'^0.3 

1 6+0* 



The table shows the confidence that the knee model can be rejected, the best fitting slope above the knee a and the best-fitting knee 
frequency i-'kncc ^^'^ 90% lower and upper confidence limits to the knee frequency (in brackets). Errors quoted for a correspond to the 
values of a below which the fit probability is reduced to less than 10% (i.e. they represent 90% confidence limits). Note that these 
confidence limits are not interpolated between sampled points in the parameter space, unlike the confidence contours plotted in Fig. |l^ 
The superscripts " and ^ mark the NGC 5506 results excluding and including the EXOSAT data respectively. 




Figure 11. The form of the high-frequency break model power 
spectrum. 



5.2.2 The high-frequency break model 

The motivation for the high-frequency break model comes 
from the power spectrum of the black hole X-ray binary 
Cyg X-1 in the low state, which shows two frequency breaks, 
as described in Section |l| If AGN have a similar power- 
spectral shape to Cyg X-1 (albeit scaled down in frequency 
by some factor), then because the power spectral slopes of 
AGN lightcurves measured at > 10"^ Hz by EXOSAT (e.g. 
Green, M'^Hardy & Lehto 1993) are significantly greater 
than 1, we may be seeing the analog of the high-frequency 
break in Cyg X-1. 

To test this possibility, we should try fitting the observed 
power spectra with a model of the form used to fit the high- 
frequency power spectrum of Cyg X-1 (e.g. Nowak 1999): 



A 



A 



otherwise 



Where A is the power-spectral amplitude at the break fre- 
quency i^bk, and Qhi and aio are the high and low- frequency 
slopes respectively, such that > aio. An example of a 
high-frequency break model with Qhi = 2 and aio = 1 is 
shown in Figure |ll|. 

Note that there is no physical basis for the sharpness of the 
break in the high-frequency break model. However, since the 
model can adequately describe the high-frequency power- 



spectral shape of Cyg X-1, it should also serve as a possi- 
ble empirical representation of the power spectra of poorer 
quality which we measure here. We do not consider the low- 
frequency break in this model in order to minimise the num- 
ber of free parameters. This approach is valid since, if the 
model is correct, low- frequency breaks will occur at least a 
decade lower in frequency than any measured high-frequency 
breaks and so will not contribute as significantly to any flat- 
tening (besides which, if additional low-frequency breaks are 
significant they will be apparent from the residuals in any 
comparison of the data with the model). 
We only consider a low frequency slope aio = 1- Clearly 
it is desirable, on grounds of computation time, to restrict 
the number of free parameters by fixing the slope below the 
break, but there are also compelling observational reasons 
why we might fix the slope to aio = 1. One particularly strik- 
ing aspect of all the Cyg X-1 power spectra is that, despite 
the variations in the position of the high and low-frequency 
breaks and the slope above the high-frequency break (e.g. as 
shown by Belloni & Hasinger 1990), the slope of the inter- 
mediate power-spectrum, between the two breaks, is always 
remarkably close to 1. Furthermore, Nowak et al. (1999) 
show that the power spectra of Cyg X-1 made from simul- 
taneous lightcurves in different energy bands show an en- 
ergy dependence above the high-frequency break (in that ahi 
decreases towards higher energies) but maintain the same 
shape (i.e. aio = 1) below the break. These results suggest 
that a power-spectral slope of 1 (or very close to 1) be- 
low the high-frequency break may, in fact, be the rule. We 
can determine if the power spectra we measure are at least 
consistent with this possibility by fitting the high-frequency 
break model (including the constant Poisson noise level, as 
before), fixing aio = 1 and stepping through the same pa- 
rameter ranges as used to fit the knee model (i.e. ahi = 1.0- 
2.4 in increments of 0.1, i^bk = 10~*-10~'' Hz in multiples 
of 2). 

The resulting best-fitting parameters are shown in Table o, 
with the results presented in the same format as for the 
knee model. The best-fitting model average power spectra 
and confidence contour plots are shown in Figure [13. The 
high-frequency break model is an acceptable description of 
the data at better than 10% confidence for all objects in 
the sample. The lower limit to break frequency is not con- 
strained in the power spectrum of NGC 5548. 



© 0000 RAS, MNRAS 000, 000-000 



16 P. Uttley, I. M. WHardy and I. E. Papadakis 




Frequency (Hz) 



l°g(^knee (H^)) 



Figure 10. Comparison of best-fitting model average power spectra (points with error bars) with observed power spectra (lines) for the 
knee model described in the text (left), and corresponding confidence contours for the parameter space searched (right). The dashed and 
solid confidence contours represent the 68% and 90% confidence limits respectively. 



Table 6. Results from fitting the broadband power spectra of four Seyfert galaxies with a high-frequency break model 



MCG-6-30-15 
NGC 5506" 
NGC 5506'' 
NGC 5548 
NGC 3516 



Rejection confidence 
33% 
10% 
3% 
27% 
39% 



See Table ^ for dcscripti( 



2.0 ±0.3 

2 1+03 
^•-^^0.7 

2 2+''-2 



i^bk/lO-*^ Hz 
51.2 (12.8-102.4) 
25.6 (0.0-102.4) 
51.2 (0.4-102.4) 
2.56 (0.0-10.24) 
2.56 (0.64-5.12) 
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Figure 12. Comparison of best-fitting model average power spectra (points with error bars) with observed power spectra (lines) for the 
high-frequency break model described in the text (left), and corresponding confidence contours for the parameter space searched (right). 
The dashed and solid confidence contours represent the 68% and 90% confidence limits respectively. 



6 DISCUSSION 

6.1 Summary of results 

We find that the power spectra of three of the four 
Seyfert galaxies in our sample (MCG-6-30-15, NGC 5506, 
NGC 3516) flatten significantly at low frequencies and that 
the power-spectral flattening can be well-fitted by either a 
knee or a high-frequency break model. Although there is no 
significant evidence for low-frequency flattening in the power 
spectrum of NGC 5548, our model fits show that we cannot 
rule out the possibility of a knee or break in the power spec- 



trum of this object also (at u < 10~® Hz ot v < 10~^ Hz for 
knee or break models respectively). 

We stress that the detection of low-frequency fiattening in 
the power spectra of AGN which we report here is com- 
pletely robust, since it is based on a rigorous Monte Carlo 
technique which we use to formally reject the simple power- 
law model for the power spectrum. On the other hand, our 
measurements of characteristic knee or break frequencies for 
the flattening are model dependent, as can be seen by the 
fact that two different models for the flattening can fit the 
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data. Clearly, the data are not yet adequate to distinguish 
between different models for the flattening. 

6.2 Comparison with naive fits to the observed 
power spectrum 

Our results confirm earlier indications of flattening in the 
power spectrum of NGC 5506 (M'^Ifardy f 989) and the ev- 
idence of flattening in the power-spectra of NGC 3516 and 
MCG-6-30-15 presented by Edelson & Nandra (1999) and 
Nowak & Chiang (2000) respectively. Chiang et al. (2000) 
also claim low-frequency flattening in the power spectrum 
of the RXTE ASM lightcurve of NGC 5548 but, as shown 
in Figure ^ ASM lightcurves of faint sources do not show 
the true source variability. All these claims of power-spectral 
flattening use a more 'naive' fit to the observed power spec- 
trum, simply fitting the assumed model directly using the 
data to determine errors at each frequency, and taking no 
account of the distorting effects of aliasing or red-noise leak. 
The fact that we confirm the power-spectral flattening re- 
ported in these previous works raises the question: is it really 
necessary to use a Monte Carlo technique to fit the observed 
power spectra? 

To answer this question, we naively fit the observed power 
spectra with the high-frequency break model, which the 
PSRESP method shows is a good fit to all the data. We bin the 
measured 'raw' power spectra (obtained from the rebinned 
and interpolated lightcurves, see Section [j.l| ) in logarithmi- 
cally spaced bins, separated by a factor of 2 in frequency 
(e.g. Edelson & Nandra 1999), using a minimum of four 
measured frequencies per bin. We determine the standard 
error from the spread in measured powers in each bin. The 
binned power spectra and best-fitting high-frequency break 
models are shown in Figure ^ and comparison of the best- 
fitting parameters with those obtained using Monte Carlo 
fits are shown in Table ^. 

The first point to note is that most of the error bars 
on the observed power spectra are much smaller than we 
would expect given the spread in power indicated by sim- 
ulated power spectra. This is because the power at each 
Fourier frequency is randomly distributed with an exponen- 
tial (xi) distribution, hence if only a few points are sampled 
the spread in points is likely to be small. A larger num- 
ber of points (> 20) must be averaged in each bin in or- 
der that error bars determined from the data are reliable 
(see Section 3.2). Because the errors are underestimated. 



the best quality broadband power spectra (for MCG-6-30- 
15 and NGC 3516, which have more extensive long-look and 
intensive monitoring lightcurves respectively) are badly fit- 
ted by the model. Better fits are obtained for the poorer 
quality data (for NGC 5506 and NGC 5548), but the 90% 
confidence errors on the model parameters are greatly un- 
derestimated, so that naive fitting implies that the flattening 
in the power spectrum of NGC 5548 is significant, whereas 
Monte Carlo fits show that it is not. 

The underestimation of errors is a major problem for naive 
model fitting, which must be taken into account when con- 
sidering claims of power-spectral flattening in the literature. 
For example, Nowak & Chiang (2000) claim a second, sig- 
nificant low-frequency break at ~ 10~^ Hz (flattening to 
zero slope) in the power spectrum of MCG-6-30-15, mea- 



sured from long-look data alone. Our simulations show that 
this additional flattening is not significant. Furthermore, our 
simulations show that the best-fitting model of Nowak & 
Chiang, when applied to the entire broadband power spec- 
trum, is ruled out at > 99% confidence: there is significant 
power at lower frequencies than 10~^ Hz. 
We now note the effects of aliasing on naive model fits to ob- 
served power spectra. To illustrate aliasing effects, we show 
the best-fitting high-frequency break models from the Monte 
Carlo fits (i.e. unmodified by aliasing or the Poisson noise 
level) in Figure |l^ as dotted lines (except for NGC 5548, 
where we show the unbroken power law model, which was 
an acceptable fit to the data). Monte Carlo fits show that 
the power spectrum of NGC 3516 is steep at high frequen- 
cies {a ~ 2). Therefore, since there is little high-frequency 
power, the amount of aliasing at lower frequencies is small 
and thus the naive fitting of a break model produces similar 
results to those obtained by Monte Carlo fits. Similarly, the 
naive fits to the power spectrum of NGC 3516 carried out 
by Edelson & Nandra (1999) yield a similar knee frequency 



to that measured using simulations (//kn 



4 X 10"'^ Hz 



versus ~ 6 x 10~^ Hz respectively). 
In contrast to the case of NGC 3516, the X-ray lightcurve of 
MCG-6-30-15 contains significant power up to high frequen- 
cies, and hence the effects of aliasing are much more signifi- 
cant, especially in the low frequency power spectrum, deter- 
mined from the long-term lightcurve which has the largest 
sampling interval. The result of this aliasing is a disconti- 
nuity in the measured broadband power spectrum from the 
medium to low-frequency parts of the power spectrum (re- 
flected in the simulated model average power spectra, see 
Figure ^ and Figure ^), so that the low- frequency power 
spectrum appears flatter and may be signficantly raised 
above the medium-frequency power spectrum. This effect is 
apparent when we compare the best-fitting model obtained 
from the naive fits with that obtained by Monte Carlo fits 
(solid and dotted lines respectively in Figure ^): the high 
frequency end of the low-frequency power spectrum is raised 
significantly above the true power level by the effects of 
aliasing, so that in the naive fits, the position of the break 
is pushed to significantly lower frequencies than those esti- 
mated by the Monte Carlo fits. The same effect might also 
cause the apparent break, where none is actually required, 
in the power spectrum of NGC 5548, since the power spec- 
trum in this case may be unbroken and have a relatively 
flat slope a = 1.6, which causes significant distortion due to 
aliasing in the high-frequency ends of both the medium and 
low-frequency power spectra. 

Note that the distorting effect of aliasing is made worse by 
the fact that the high-frequency ends of the power spectra 
are also the most well-defined, as they are made by aver- 
aging over a large number of frequencies, so that system- 
atic errors due to aliasing are more pronounced than they 
would be if the low-frequency ends of the power spectra were 
well sampled (since for red-noise power spectra, the aliased 
power at frequencies much lower than the Nyquist frequency 
is small compared to the underlying power). Because of this 
effect, particular caution should be applied to claims of low- 
frequency flattening based on naive fits to the power spectra 
of AGN which are highly variable on short time-scales (e.g. 
Akn 564, Pounds et al. 2001). 

Under what circumstances might naive fitting produce re- 
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Figure 13. Best-fitting high-frequency break models (soUd lines) fitted naively to binned observed power spectra (points with error bars). 
Error bars are standard errors determined from the data. To illustrate aliasing effects we also show the best-fitting underlying power 
spectral models obtained from the Monte Carlo fits (dashed lines), using the broken power law model for all sources excpet NGC 5548, 
where we show the simple power law (see text for discussion). 



Table 7. Comparison of naively fitted break model parameters with Monte Carlo results 

Naive results Monte Carlo results 

X^/d.o.f. a i^bk/lO"*^ Hz a i^bk/lO"*^ Hz 

MCG-6-30-15 59.6/22 1.68 5.43 2.0 ± 0.3 51.2 (12.8-102.4) 

NGC 5506 25.2/20 1.95±j;-2| 20.0 (9.2-44.2) 2.4±j; g 51.2 (0.4-102.4) 

NGC 5548 17.2/11 1.9±g-^** 0.87 (0.18-1.9) 2.4±% 2.56 (0.0-10.24) 

NGC 3516 60.16/22 2.06 1.9 2.2±0-2 2.56 (0.64-5.12) 



The table shows the x (and degrees of freedom), slope above the break, a and break frequency, i/bk obt a.iried 



and for comparison the model parameters estimated from the Monte Carlo fits using PSRESP (see Section 5.2.2). Errors are 90% 
confidence, corresponding to Ax^ = 4.61 in the naive fits. Errors are not quoted where the fit is very bad. 



from naive model fits. 



liable results (and associated uncertainties)? Clearly, the 
eflFects of aliasing can be mitigated by binning up the 
lightcurves to a longer sampling interval, thus smoothing 
out the variability on time-scales less than the sampling in- 
terval used to make the power spectrum. Unfortunately, the 
effect of this binning is to reduce the Nyquist frequency and 
hence the frequency range sampled by the power spectra. 
Furthermore, since reliable errors may only be estimated by 
averaging over < 20 points per frequency bin, the lowest fre- 
quency sampled in our observed power spectra would have 
to be increased by a factor ~ 10, in order that sufficient 
low-frequency power-spectral points are averaged over. The 
net effect of binning up the lightcurves and power spectra is 
to reduce the frequency range covered by the individual low 
and medium-frequency power spectra by ~ 2 decades (i.e. 
to virtually a single frequency bin in each), while the low- 
frequency end of the high-frequency power spectrum (which 
does not suffer so much from aliasing) is shifted up by a 
decade. 

Therefore, we conclude this discussion by noting that the 
naive method may be used to fit power spectra over rela- 
tively narrow ranges which are well sampled (i.e. the highest 
frequencies in AGN long- look lightcurves), but the Monte 



Carlo method of power-spectral fitting, of which PSRESP is 
an example, is essential in order to place reliable constraints 
on the underlying broadband power-spectral shape of AGN 
lightcurves obtained to date. 

6.3 Estimating black hole masses 

Having measured characteristic knee or break frequencies 
under the assumption of different underlying power-spectral 
models, we can consider the implications of the measured 
characteristic frequencies for the black hole mass. Through- 
out our discussion, we make the specific assumption that 
all the objects in our sample have the same power-spectral 
form (knee or break) so that the characteristic frequencies 
are directly comparable. 

What is most striking about the characteristic frequencies 
measured for MCG-6-30-15 and NGC 3516 is that they are 
significantly different, in either model, at a level of better 
than 99% confidence, i.e. the 90% confidence limits of the 
break frequencies do not overlap. This is unexpected, be- 
cause both objects have similar 2-10 keV X-ray luminosities 
of 1/2-10 ~ 1.5 X 10^^ erg s~^, so the fundamental parame- 
ter driving the position of the characteristic frequency is not 
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Table 8. Black hole mass and bolometric luminosity estimates 

knee model high-frequency break model 

best estimate upper mass limit best estimate upper mass limit 

M/106Mq Lboi/iEdd M/106Mq Ltoi/LjsM M/W^Mq Lboi/^Edd M/W^^Mq Lboi/^-Edd 

MCG-6-30-15 0.25 2.4 1.6 0.4 0.6 1 5 0.14 

NGC 5506 0.5 6 25 0.1 0.6 5 150 0.02 

NGC 5548 65 0.2 NA NA 13 0.8 NA NA 

NGC 3516 2 1.5 13 0.25 13 0.25 94 0.03 

See text for details. Mass estimates for NGC 5506 use the model fit results for the RXTE plus EXOSAT data. No upper limits can be 
set on the black hole mass of NGC 5548. 



strictly related to the luminosity. One intriguing possibility 
is that although both objects have a similar luminosity, they 
may have different black hole masses. It is quite possible that 
the break or knee time-scales scale linearly with black hole 
mass, as would be expected if they correspond to a charac- 
teristic time-scale in the accretion disk or the characteristic 
size scale of the system. If this is the case, then MCG-6-30-15 
must have a significantly smaller black hole than NGC 3516 
and so must be accreting at a much higher fraction of its 
Eddington limit. 

Let us assume that the characteristic frequency scales lin- 
early with the inverse of black hole mass, and that the mass 
of the black hole in Cyg X-1 is 10 Mq (Nowak et al. 1999). 
We will scale the best-fitting frequency measured by the knee 
model with the mean frequency of the low-frequency break 
in the low-state power spectrum of Cyg X-1, which varies 
between 0.04-0.4 Hz but has an average value of 0.13 Hz 
(taken from the 28 separate measurements presented in Bel- 
loni & Hasinger 1990). Similarly, we scale the best-fitting 
frequency measured in the high-frequency break model with 
the high-frequency break in the low-state power spectrum 
of Cyg X-1, with mean value 3.3 Hz, varying between 1- 
6 Hz. The resulting best estimates of black hole mass are 
shown in Table ^, together with conservative upper limits to 
the black hole mass, obtained by scaling the 90% confidence 
lower limits to knee/break frequency with the respective up- 
per values of the observed knee/break frequency ranges in 
Cyg X-1 (i.e. 0.4 Hz and 6 Hz for knee and break frequencies 
respectively). We also show the corresponding bolometric 
luminosities expressed as a fraction of the Eddington lumi- 
nosity for a black hole with the estimated mass, LEdd — 
1.3 X IO^^Mbh erg s"^ For MCG-6-30-15, we assume the 
bolometric luminosity of 8 x 10^^ erg s~^, estimated from a 
multiwavelength study by Reynolds et al. (1997). For the re- 
maining AGN, bolometric luminosities are derived from the 
typical 2-10 keV X-ray luminosities (1.5 x 10*^ ergs s~^ for 
NGC 5506 and NGC 3516, 5 x lO"*^ ergs s"^ for NGC 5548) 
by applying the mean bolometric correction of Padovani & 
Rafanelli (1988), Lboi — 271/2-10- Note that scatter in this 
bolometric correction for individual sources, estimated from 
the measurements of Padovani & Rafanelli (1988) is possibly 
up to a factor ~ 3 in either direction. The bolometric cor- 
rection for MCG-6-30-15 is a factor 5 lower than the mean 
value, but Reynolds et al. (1997) note that their estimate 
is a lower limit, based on the assumption of minimal red- 
dening in this source, so that the bolometric luminosity in 
MCG-6-30-15 may be substantially larger than shown here. 
Table ^ shows that, for the objects with the best-constrained 
characteristic frequencies in their power spectra (MCG-6-30- 
15 and NGC 3516), the knee model predicts a lower black 



hole mass than the high-frequency break model. Reverbera- 
tion mapping estimates of the mass of NGC 3516 suggest a 
black hole mass of ~ 2 x 10^ Mq (Wanders & Home 1994), 
compatible with the best estimate of mass from the high- 
frequency break model, but not inconsistent with the upper 
limit derived from the knee model. Given the uncertainty in 
the true bolometric luminosity, and the conservative upper 
mass limits quoted here, the masses predicted by the knee 
model are consistent with sub-Eddington accretion rates. 
However, the conservative 40% Eddington lower limit to ac- 
cretion rate estimated from the knee model fits to the power 
spectrum of MCG-6-30-15 (where the bolometric luminosity 
quoted represents a firm lower bound) is in contrast to that 
observed in Cyg X-1 where the transition between the low 
and high states occurs at only a few per cent of the Edding- 
ton luminosity (Di Salvo et al. 2001, Zhang et al. 1997). Since 
the knee model represents the low-frequency break seen in 
the power spectrum of Cyg X-1 in the low state, it would 
seem that the evidence from accretion rate is inconsistent 
with the knee model, at least in this source. Interestingly, 
even the high-frequency break model predicts an upper limit 
to black hole mass which is only barely consistent with the 
low-state accretion rates seen in Cyg X-1. An intriguing 
possibility is that MCG-6-30-15 is in a state analogous to 
the high state of Cyg X-1, which has a characteristic steep 
(q ~ 2) high-frequency power spectrum breaking at 10 Hz 
to Q ~ 1, with no low- frequency breaks down to frequencies 
as low as 10~^ Hz (Cui et al. 1997). Therefore, high state 
power spectra similar in shape to that of Cyg X-1 would 
also be compatible with the high-frequency break model we 
use here. Scaling by the high state frequency-break yields 
an estimated black hole mass of 2 x 10® Mq (upper limit 
8 X lO" M0), implying an accretion rate of ~ 30% Ed- 
dington (which we might expect from a high-state source). 
Although the black hole mass in MCG-6-30-15 has not yet 
been measured by reverberation mapping, Reynolds (2000) 
has pointed out that the low luminosity (absolute B magni- 
tude ~ —19) of the SO host galaxy of this AGN corresponds 
to a black hole mass of ~ 10'' M0 if the black hole mass- 
bulge luminosity relation is the same as that determined by 
Magorrian et al. (1998) for a sample of normal elliptical and 
SO galaxies. We note here that the black hole masses of all 
the Seyfert galaxies which have been reverberation mapped 
show smaller black hole masses by a factor of ~ 10 than 
the Magorrian relation suggests (Wandel 1999). This dis- 
crepancy would imply a black hole mass for MCG-6-30-15 
which is significantly lower than 10^ Mq, consistent with 
the low mass suggested by the high state interpretation of 
the power spectrum. We further note that the lower, rever- 
beration mapping estimates of black hole mass in AGN also 
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follow the recently discovered correlation of black hole mass 
with bulge velocity dispersion, which has been found to be a 
better predictor of black hole mass than the Magorrian rela- 
tion (Gebhardt et al. 2000, Ferrarese et al. 2001). Therefore 
it seems plausible that the black hole in MCG-6-30-15 has a 
mass of order 10^ M©. Reverberation mapping, or another 
independent mass estimation technique will be required to 
confirm this. 

The power-spectral shape of NGC 5506 is also consistent 
with a high state interpretation of this source, although the 
uncertainties on the break frequency are sufficiently large 
that a low-state interpretation would suffice to fit the data 
self-consistently. We note here however that the low luminos- 
ity of the NGC 5506 galaxy (absolute B magnitude ~ —20), 
suggests that the central black hole mass is relatively low. 
Although reverberation mapping of NGC 5506 is not possi- 
ble because it is a Seyfert 2, it does have a water megamaser 
(Braatz, Wilson & Henkel 1994) which might permit a good 
estimate of its central black hole mass for comparison with 
the mass estimated from its power spectrum. 
The high state AGN interpretation for either MCG-6-30-15 
or NGC 5506 is highly speculative at this stage. One signif- 
icant problem is that the energy spectra of neither source 
shows the steep power-law (F ^ 2.6) and strong disk black- 
body components seen in the spectrum of Cyg X-1 in the 
high state (Cui et al. 1998), although this difference might 
be a result of the significantly lower disk temperature ex- 
pected in AGN compared to BHXRBs. Greater support for 
the high state interpretation in either MCG-6-30-15 and 
NGC 5506 could be provided by the detection of significant 
power down to very low frequencies, implying that there is 
no low-frequency break in the power spectrum on frequen- 
cies < O.Olfbk. Continued long-term monitoring will help to 
resolve this issue. 

The frequency break in NGC 5548 is not signficantly de- 
tected. The black hole mass of NGC 5548 measured from 
reverberation mapping is ~ 10* Mq (Wandel, Peterson & 
Malkan 1999), implying that the accretion rate is relatively 
low (< 10% Eddington), so applying the low state interpre- 
tation we can estimate that the high-frequency break occurs 
at ~ 2 X lO"'^ Hz, and may become detectable with further 
monitoring observations. 



6.4 The AGN-BHXRB connection and physical 
implications 

The main result of this work has been to show that the power 
spectra of AGN fiatten towards low frequencies. Whether 
this flattening has the same form as that in BHXRBs is 
still not certain. The data do not yet allow us to categori- 
cally rule out either knee or high-frequency break models, let 
alone more complex models where the power spectra might 
have multiple breaks or flatten gradually. However, it is en- 
couraging that a simple high-frequency break, similar to that 
seen in the classic BHXRB Cyg X-1, can fit the data well, 
and yield plausible black hole masses if we assume that the 
break frequency scales inversely with black hole mass. 
A useful way of comparing power spectra is to plot 
frequency X power, rather than power, as a function of fre- 
quency. The vPiy) plot produced in this way is analogous to 
the vFiy) method of displaying energy spectra, in that the 



peak in the vPiy) plot shows which time-scales most of the 
variability occurs on, as well as the magnitude of variability 
on those time-scales. In Figure ^ we plot the h'P{u) power 
spectra for the best-fitting knee and high-frequency break 
models fitted to the observed power spec tra of MCG-6-30- 
15, NGC 5506 and NGC 3516 (see Section ^.2. 2|). The uP{u) 



power spectra of NGC 5548 are also included in the figure, 
except that since the knee or break is not defined for this ob- 
ject, we have assumed knee and break frequencies of 10~* Hz 
and 2 x 10~^ Hz respectively, corresponding to the black hole 
mass of 10* M© estimated from reverberation mapping. For 
comparison, vP{u) power spectra of Cyg X-1 are included, 
corresponding to typical power-spectral parameters in the 
low and high states (Nowak et al. 1999, Cui et al. 1997). 
Note that the flat peaks in the high-frequency break vP{y) 
power spectra of AGN and the typical power spectra of 
Cyg X-1 correspond to the a = 1 part of the power spec- 
trum, where there is equal integrated power per decade of 
frequency. The low-frequency drop-off in the vP{v) power in 
the low state of Cyg X-1 corresponds to the low- frequency 
break or 'knee' in the power spectrum, which we have not 
yet detected in the AGN power spectra if we are seeing the 
high-frequency break, because our low-frequency data is not 
yet adequate to detect it. Therefore, until clear evidence for 
low- frequency flattening to a slope a = 1 is found, we can- 
not rule out the possibility that all our sources have power 
spectra similar to the high state power spectrum in Cyg X- 
1. 

Note that the knee model yPiy) power spectra do not 
show flat peaks, because the high-frequency slopes fitted 
are all significantly steeper than 1. If we were measuring the 
low-frequency break alone (i.e. the frequency of the high- 
frequency break is too high to detect), we would expect the 
knee model to fit high-frequency slopes of a ~ 1. This re- 
sult provides additional evidence that we are indeed measur- 
ing the high-frequency breaks in these sources, if the power 
spectra of AGN have the same shape as those of BHXRBs. 
However, we stress that we cannot yet rule out the possi- 
bility that we are measuring more complex configurations 
including both high and low-frequency breaks. 
Fig. ^ shows clearly that the vP{v) power spectra of AGN 
are similar to those in Cyg X-1, in that they have similar 
peak powers, but are shifted at least 5 decades down the 
frequency axis. The fact that the peak powers are similar 
implies that the number of varying regions and the general 
pattern of variability is the same in AGN and BHXRBs, 
while the luminosity of the varying regions and their vari- 
ability time-scales are scaled by some factor which could be 
the black hole mass. 



7 CONCLUSIONS 

We have presented long-term RXTE monitoring lightcurves 
for 4 Seyfert galaxies, MCG-6-30-15, NGC 5506, NGC 5548 
and NGC 3516, and measured their broadband power spec- 
tra to determine if they flatten towards low frequencies, like 
those of BHXRBs. The interpretation of power spectra mea- 
sured from discretely (and sometimes unevenly) sampled 
lightcurves is complicated by the distorting effects of red- 
noise leak and aliasing. A further serious problem is that 
reliable errors in the power in each frequency bin cannot be 
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Figure 14. Comparison of best-fitting knee (left panel) and high-frequency break (right panel) model vP{y) power spectra of MCG-6- 
30-15, NGC 5506, NGC 3516, NGC 5548 with power spectra of Cyg X-1 in the low and high states. 



defined from the data, due to the small number (< 20) of 
power-spectral measurements in each frequency bin. 
To overcome those difficulties, we have built on the response 
method of Done et al. (1992), to develop a reliable Monte 
Carlo method for testing models for the underlying power- 
spectral shape of discretely, unevenly sampled lightcurves. 
Our method, psresp, takes proper account of the effects of 
aliasing and red-noise leak and crucially, uses the distribu- 
tion of simulated power spectra to define reliable confidence 
limits on our model fits. 

We have used psresp to test simple models for the power- 
spectral shape of the active galaxies in our sample. Our main 
results are: 

1. A single power-law model for the power spectra, with 
no low-frequency flattening, is rejected at better than 95% 
confidence for MCG-6-30-15, NGC 5506 and NGC 3516. The 
power spectrum of NGC 5548 is consistent with a single 
power-law. 

2. Both the knee model (flattening to slope a = 0) and 
high-frequency break model (flattening to a = 1) pro- 
vide good fits to the power spectra of all four sources. 
Knee/break frequencies are well constrained for all sources 
except NGC 5548, for which we can define upper limits only. 

3. The characteristic knee/break frequency measured for 
MCG-6-30-15 is significantly higher (99% confidence) than 
the corresponding frequency in NGC 3516, even though 
both sources have similar X-ray luminosities (~ 1.5 x 
lO"*^ erg s"^). 

4. If the knees or breaks correspond to those seen in the 
low-state power spectrum of Cyg X-1, and the characteristic 
knee/break time-scales scale linearly with black hole mass, 
the black hole masses estimated from the knee model are 
lower than those estimated by the break model, but remain 
consistent with sub-Eddington accretion rates. In the case 
of MCG-6-30-15, the conservative lower limit on accretion 
rate estimated from the knee model is an order of magnitude 
higher than that seen in the low state of Cyg X-1, favour- 
ing the high-frequency break model and further suggesting 
that the break seen in the power spectrum of MCG-6-30-15 



(and possibly NGC 5506) is analogous to the high-frequency 
break seen in the high state power spectrum of Cyg X-1. 

5. The pPiy) power spectra of the Scyfcrt galaxies stud- 
ied here are similar in amplitude to the uP{v) power spectra 
of Cyg X-1 in the high and low states. 

Conclusions 4 & 5 imply that the power spectra of AGN 
are consistent with being identical in shape and fractional 
RMS amplitude (integrated over the whole power spectrum) 
to those of BHXRBs, with characteristic time-scales scaling 
linearly with black hole mass. Arguments based on accre- 
tion rate seem to favour the high-frequency break model 
over the knee model in at least one source (MCG-6-30-15), 
although this evidence is circumstantial. Further monitoring 
observations (which are currently underway) are needed to 
reject either model on the basis of power-spectral measure- 
ments alone, and to confirm the interesting possibility that 
two of the objects in our sample are analogous to high-state 
BHXRBs. 
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